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Abstract 

Effective Field Theory(EFT) is, the unique, model independent and systematic low-energy version of QCD 
. .for processes involving momenta below the pion mass. A low-energy photo-nuclear observable in three-body 
systems, photon polarization parameter at thermal neutron energies is calculated by using pionless EFT up 
to next-to-next to leading order (N 2 LO). In order to make a comparative study of this model, we compared 
our results for photon polarization parameter with the realistic Argonne vi$ two-nucleon and Urbana IX or 
£>5 Tucson-Melbourne three-nucleon interactions. Three-body currents give small but significant contributions to 
<^ some of the observables in the neutron-deuteron radiative capture cross section at thermal neutron energies. 

In this formalism the three-nucleon forces are needed up to N 2 LO for cut-off independent results. Our result 
CN converges order by order in low energy expansion and also cut-off independent at this order. 
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I. INTRODUCTION 



The study of the three-body nuclear physics involving nucleon-deuteron photodisintegration of 3 He 
and 3 H as well as the time reversed nucleon radiative capture by deuteron has been investigated 
in theoretical and experimental works over the past decays. The photon polarization in the reaction 
neutron-deuteron capture reaction has been measured with polarized thermal neutrons by Konijnenberg 
et al. [if. In addition the search for three- nucleon force effects in the electromagnetically induced 
process, come more and more into the focus in the recent years. In principle, using the continuity 
equation for three-nucleon forces lead to three-nucleon currents. It is a quantitative question based 
on current choices of nuclear force models to reveal signatures by switching on and off three-nucleon 
forces. 

Several groups are studying electromagnetic processes in the three-nucleon system. In Refs. 0, [sj , the 
nucleons are taken as interacting via two- and three-nucleon potentials and the electromagnetic currents 
are then constructed using the exchange scheme to satisfy the current conservation relation(CCR), but 
only with a part of the interaction. In another work, the meson exchange currents are taken into 
account using the Siegert's theorem [ij without three-body currents. Three-body currents are added 
using a nuclear model which allows for the excitation of a nucleon to a A isobar. The A excitation 
yields also effective three-body forces and three-body currents. However, this current model does not 
satisfy exactly the CCR with the adopted Hamiltonian [51] . In alternative descriptions of the three- 
nucleon electromagnetic processes for very low energy pd and nd radiative capture, in model dependent 
theory, a variety of electromagnetic observables involving the two- and three-nucleon forces have been 
extensively studied in the past by several research groups (for a review, see Ref. 0). Recently, Viviani 
et al. have investigated the nd radiative capture reactions below deuteron breakup threshold [7]. Their 
work shows sensitivity to short-range physics namely, details of including the physics of the Delta and 
pion-exchange currents. They obtained the cross section from Argonne vn two-nucleon and Urbana 
VIII three-nucleon interactions (AV14/UVIII), also from Argonne Vx$ two-nucleon and Urbana IX three- 
nucleon interactions(AV18/UIX), including A admixtures. They found Cross section of 0.600 (mb) and 
0.578 (mb) which are above the experimental values by 18% and 14%, respectively. It is worthy to 
mention that the explicit inclusion of A-isobar degrees of freedom in the nuclear wave function improve 
the agreement with the experimental data better than those obtained using the perturbation theory, 
Apy. This indicate that their results for very-low energy observables are sensitive to the details of 
the short-range part of the interaction. Recent calculations using gauge-invariant currents reduced the 
spread Q|, however including three-body currents results 0.558 mb, which still is above the data by 
10%. Model-dependent currents associated with A(1232) were identified as a source of the discrepancy. 
Thus, the question that remains is, how such details of short-range physics can so severely influence a 
very long range reaction with maximal energies of less than 10 MeV. 

Recently developed pionless EFT is particularly suited for high order precision calculation. The 
so-called pionless EFT in nuclear physics aspires a systematic classification of all forces. At its heart 
lies the tenet that physics at those very low energies can be described by point-like interactions between 
nucleons only. In this approach, all particles but the nucleon themselves are considered high energy 
degrees of freedom and are consequently "integrated out". The resulting EFT is considerably simpler 
than potential models or the "pionful" version of nuclear EFT (in which pions are kept as explicit 
degrees of freedom), but its range of validity is reduced to typical momenta below the pion mass. 
There are many processes situated at thermal energies which are both interesting in their own right 
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and important for astrophysical applications. Recently we have calculated the cross section of radiative 
capture process nd — > 3 if 7 by using pionless EFT [3, llOj. No new three-nucleon forces are needed 



up to N 2 LO in order to achieve cut-off independent results, in addition to those fixed by the triton 
binding energy and nd scattering length in the triton channel. The cross-section is determined to be 
a tot = [0.503 ± 0.003]m6. 

The present study investigates a low-energy photo-nuclear observable in three-body systems namely, 
photon polarization parameter at thermal neutron energies, using pionless EFT up to N 2 LO. The 
emphasis is on constructing three-body currents with model independent theory corresponding to three- 
nucleon interactions and comparison of the our model's result with those of other model dependent 
theory. 

The paper is organized as follows: In Section [Til we briefly review theoretical framework including 
the Faddeev integral equation, three-body forces and cut-offs dependence for calculating of neutron- 
deuteron radiative capture at thermal energies. Then we calculate the photon polarization parameter at 
thermal energies in Section lTTTI and compare our results with the corresponding experimental. Section HVl 
is devoted to comparison of our result with the data and the results of other theoretical models. 
Summary and conclusions are given in Section M 



II. NEUTRON-DEUTERON RADIATIVE CAPTURE AT THERMAL ENERGIES 

At thermal energies the nd capture reaction proceeds through S'-wave capture predominantly via 
magnetic dipole transitions from the initial doublet J=l/2 and quartet J=3/2 nd scattering states. 
In addition, there is a small contribution due to an electric quadrupole transition from the initial 
quartet state. Consequently, 2 Si/ 2 describes the preferred mode for nd — > 3 Hj and pd — > 3 He r y. The 



three-nucleon Lagrangian is well-known and will not be discussed here [bH TEI 



As long-distance phenomena must be insensitive to details of the short range physics (and in par- 



ticular of the regulator chosen), Bedaque et al. [14], [15( showed that the system must be stabilized by 
a three-body forces 



H(E; A, = £ f>.(A) (^) " = ^ + ^ (ME + + . . . . (!) 

n=0 ^ ' 



which absorbs all dependence on the cut-off as A — > oo. Eq.([T]) is analytical in E and can be obtained 
from a three-body Lagrangian, employing a three-nucleon auxiliary field analogous to the treatment 
of the two-nucleon channels [l4j. Contrary to the terms without derivatives, the term involves three- 
body forces (second term) contains two derivatives. The derivation of the integral equation describing 



neutron-deuteron scattering has been discussed before [15] . We present here only the result, including 
the new term generated by the second term in Eq.flT]). The resulting amplitudes is a mixture of t s 
describes the d t + N — > d s + N process, and t t describes the d t + N — > d t + N process: 
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t s (p,k) = l l [3)C( P ,k) + 2H(E,A)} + ±- J dqq 2 [V s (q)[]Cip,q) + 2n(E,A)]t s (q) 

o 

+V t (q) [3JC(p,q) + 2H{E,A))t t (q)} 

A 

t t (p, k) = \ [JC(p, k) + 2H(E, A)] + i- J dq q 2 [V t (q) [JC(p, q) + 2H(E, A)] f t (g) 

o 

+V s (q) [3/C(p, q) + 2H(E, A)] , (2) 

2 

where T> Sjt (q) = T> s ^{E—^r, q) are the propagators of deuteron. For the spin-triplet S-wave channel, one 
replaces the two boson binding momentum 7 and effective range p by the deuteron binding momentum 
7i = 45.7025 MeV and effective range p t = 1 .764 fm. Because there is no real bound state in the spin 
singlet channel of the two-nucleon system, it is better to determine the free parameters by the scattering 
length a s = l/jg = —23.714 fm and the effective range r s = 2.73 fm at zero momentum QEH. The 
neutron-deuteron J = 1/2 phase shifts 5 is determined by the on-shell amplitude t t (k,k), multiplied 
by the wave function renormalisation 

3-7T 1 

m = Ztt{k ,k) = - I ^— k . (3) 

The spine structure of the matrix elements for neutron radiative capture by deuteron is complicated, 
however in very low energy for this reaction we can introduced three multipole transition that is allowed 
by p-parity and angular momentum conservation i.e. I p = | — ► Mi and I p = | — > Mi,E 2 . The 
parameterization of the corresponding contribution to the matrix elements and the Ai\ amplitude are 
from the magnetic moments of the nucleon and dibriyon. These are well-known and will not be given 
here @,[Ioj. 

The radiative capture cross section nd — > 3 Hj at very low energy is given by 0] , 

Hr^Eii^'il (4) 

y v rel zm n . ls j 

where 

xf SJ = —V^xt SJ (5) 
PPn 

where \ is either the magnetic or electrical moment and p^ is nuclear magneton and p is momentum 
of the incident neutron in the center of mass. The contribution of the electric transition E[ J to total 
cross section at energies less than 60 KeV is insignificant. Therefore, the electric quadrupole transition 
^0(3/2) (3/2) f rom j n jtj a } quartet state will not be considered at thermal energies. 

We now turn to the Faddeev integral equation used in the magnetic moment calculation and also 
the interaction kernel included in this integral equation. Fig. [1] represents the contribution diagrams 
for adding photon- interaction to the Faddeev equation (j2J). In these diagrams photon is minimally 
coupled to nucleons in three-body systems. The diagrams for adding photon-interaction to the Faddeev 
equation up to N 2 LO are depicted in Fig. [2j Photon is coupled to two-body system via Li vertices. 
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FIG. 1: (Color online) Some diagrams for adding photon-interaction to the Faddeev equation up to N 2 LO. 
Thick solid line is propagator of the two intermediate auxiliary fields D s and Dt, denoted by D; /C: propagator 
of the exchanged nucleon. Photon is minimally coupled to nucleons in three-body systems. Wavy line shows 
photon and small circles show magnetic photon interaction. 



Two-Body: 




FIG. 2: (Color online) Some diagrams for adding photon-interaction to the Faddeev equation up to N 2 LO. 
Photon is minimally coupled to two-body system and three-body vertices in three-body systems. For L\ 
vertices, see Ref [la]; i^three- body force, see eq.(l). Remaining notation as in Fig. [TJ 



The coefficient L% is fixed at its leading non- vanishing order by the thermal cross section [l6[. We 
have other possible diagram that can be considered for our calculation for inclusion of photon to the 
three-body vertices TC. This diagram is shown in Fig. [2J 

All corrections contribute to observables typically as Q n = Oj^j compared to the LO result and 

that low-energy observables must be independent of an arbitrary regulator A up to the order of the 
expansion. In other words, the physical scattering amplitude must be dominated by integrations over 
off-shell momenta q in the region in which the EFT is applicable, q < A^. Typical low momentum scales 
p typ in the three-body system are the binding momenta of the two-nucleon real and virtual bound states, 
7 S ~ —8.0 MeV, 7j m 45 MeV and the scattering momentum k. In addition, the three-body forces are 
determined in part by the typical three-nucleon bound state momentum 7^ ~ y/MBd ~ 90 MeV, Bd 
is the triton binding energy. The breakdown scale A^ « m n of the theory is the scale at which higher 
order corrections become comparable in size. One can therefore estimate sensitivity to short- distance 
physics, and hence provide a reasonable error analysis, by employing a momentum cut-off A in the 
solution of the Faddeev equation and varying it between the breakdown-scale A^ to 00. If observables 
change over this range by "considerably" more than Q n+1 , a counter-term of order Q n should be added. 
This method is frequently used to check the power counting and systematic errors in pionless EFT with 
three nucleons, see e.g. most recently [15l. A similar argument was also developed in the context of 
the EFT "with pions" in nuclear physics [171, |l_8 . 
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TABLE I: Comparison between different theoretical and experimental results for Neutron radiative capture 
by deuteron at zero energy (0.0253 ev). Last row shows N 2 LO order pionless EFT result. 



Experiment 


Year 


Total cross section(mb) 


Jurney et.al. [Ill 


1963 


0.60 ± 0.05 


Merritt et.al. [12] 


1968 


0.521 ± 0.009 


Jurney et.al. [13] 


1982 


0.508 ± 0.015 


Theory 






AV14/VIII(IA+MI+MD+A) [7] 


1996 


0.600 


AV18/IX(IA+MI+MD+A) [7] 


1996 


0.578 


AV18/IX (gauge inv.) [8] 


2005 


0.523 


AV18/IX (gauge inv. + 3N-current) [8] 


2005 


0.556 


EFT(N 2 LO)+3N-forces [10] 


2006 


0.503 ± 0.003 



AV14 + UVIII X X X 

no A pert. A full A 

avis + u ix -m 

exp 

EFT(X) n ° nsw 3bf 

i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — 

0.55 0.6 0.65 Gfmbf 



FIG. 3: (Color online) Comparison between calculated cross section of neutron radiative capture by deuteron 
by different theoretical models, pionless EFT and experimental data. 



III. PHOTON POLARIZATION PARAMETER AT THERMAL ENERGIES 

If the process is dominated by S-wave capture, as in the case for the neutron-deuteron radiative 
capture reaction at thermal neutron energies, the observable for circular polarization Pr{0) is simply 
given by: 

P r (6) = R c P N -q (6) 

where P/v is the polarizations of the spin-1/2 nucleon and R c is the polarization parameter (for more 
detail see Q). This polarization parameter depends on the relative sign between the amplitudes 1/2 
and 3/2 channels. Numerically, R c lies in the region —1/2 < R c < 1. Its experimental value is R c = 
-0.42 Q. 
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IV. RESULTS 



We numerically solved the Faddeev integral equation up to N 2 LO. We used he = 197.327 MeVfm, 
and mass of M — 938.918 MeV for nucleon. A deuteron binding energy (momentum) of B = 2.225 MeV 
(jd = 45.7066 MeV) is used for the nucleon-nucleon triplet channel. A residue of Z d = 1.690(3) is used 
for the NN singlet channel. The 1 So scattering length is chosen to be a s = —23.714 fm. After fixing 
the leading non- vanishing order in the thermal cross section L\ is found to be -4.5 fm. 

As in Ref. we determined which three-body forces are required at any given order, and how 
they depend on the cut-off. Low energy observables must be insensitive to the cut-off, namely to any 
details of short range physics in the region above the break down scale of the pionless EFT ( which set 
approximately by the pion mass ). 

The results for the thermal energy cross section and photon polarization parameter are presented in 
table I and II, along with the experimental data [j], 13]. Table I compares the nd ^ 3 Hj cross section 



at zero energy (0.0253 eV) for various experimental and theoretical works. The corresponding values 
for the cross section from the pionless EFT evaluation up to N2LO is shown in the last row. The EFT 
results for this cross section are presented only up to three significant digits. 

Recently in a model dependent two-body current calculation, the total cross section for nd — > 3 is 
obtained as cr-p = 0.523 mb This value can be compared with the corresponding result cr-r = 0.558 
mb obtained in Ref. [7J]. The later work used the present MI two-body current operators therefore 
leads to an estimate closer to the experimental data cr-r = 0.508 ± 015 mb [Hjj]. However, the addition 
of the three-body currents, which give a rather sizable contribution as can be seen from the row 
labeled "full-new" in table I, brings the total cross section to cr-r = 0.556 mb. Table I shows also 
EFT result of the cross section for this reaction up to N 2 LO order. Three-nucleon forces are needed 
up to N 2 LO order for cut-off independent results. Hence the cross-section is in total determined as 
a tot = ]0.485(LO) + 0.011 (NLO) + 0.007(A 2 LO)] = [0.503 ± 0.003]m6. The theoretical accuracy may 
for example be estimated conservatively by Q ~ — ~ | of the difference between the NLO and N 2 LO 
results. 

Table II shows Comparison between the results of different models dependent, model independent 
EFT and experiment for the photon polarization parameter. The photon polarization parameter is 
sensitive to two-body currents ( for its definition in terms of RMEs, see Ref. 0]). We compare our 
prediction for the photon polarization parameter with the theoretical and the experimental results 
of Ref. [l|, 0, [s| in this table. The magnetic Mi-transition gives the dominant contribution for our 
calculation. 

In Fig. [3], we compare our results with those obtained in Refs. 0, El] • As can be seen by inspecting 
Fig. [3X the pionless EFT calculations is converges order by order in low energy expansion and also 
cut-off independent at this order (see [13]). There are no new three-nucleon forces besides those already 
fixed in nd scattering at the same order. The contribution from the photon coupling to a three-nucleon 
force is negligible in our calculation. Our calculation has a systematic error which is now smaller than 
the experimental error-bar. 

V. SUMMERY AND CONCLUSION 

We applied pionless EFT to find numerical results for the photon polarization parameter Rq- At 
very low energies, the interactions between nucleons can be described only by point-like interactions. 
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TABLE II: Comparison between different theoretical results for photon polarization parameter Rc of the 
reaction neutron-deuteron radiative capture at thermal energies. The last line quotes deviation between data 
and theory, if it is larger than the theoretical or experimental uncertainty. 



Type of theory and experiment 


Rc 


Overestimation the experimental 


AV14/VIII(IA+MI+MD+A) [7] 


-0.420 


below 1% 


AV18/IX(IA+MI+MD+A) [7J 


-0.469 


12% 


AV18/IX (gauge inv.+ 2N-current) [8] 


-0.429 


2% 


AV18/IX (gauge inv. + 3N-current) [8] 


-0.476 


13% 


EFT(LO) 


-0.387 


8% 


EFT(NLO) 


-0.403 


4% 


EFT(N 2 LO) 


-0.412 


2% 


Experiment [1] 


-0.42 ±0.03 





One cannot identify pions as the lightest exchanged particles between nucleons as long as the typical 
external momentum p typ in a reaction is below the pion mass m n . That is because the Compton 
wavelengths are not small enough to resolve the nuclear forces as originating in part from one pion 
exchange. Then all particles but nucleons are integrated out. One can identify a small, dimensionless 
parameter Q = <C 1, where Kf ~ m n is the typical momentum scale at which the one pion exchange 

is resolved and pionless EFT must break down. Incident thermal neutron energies have been considered 
for this capture process. 

The photon polarization parameter R c of the reaction neutron-deuteron radiative capture nd — > 
at thermal energies was calculated in pionless EFT. This model independent and systematic low energy 
version of QCD is suited for processes involving momenta below the pion mass. At these energies 
our calculation is dominated by S'-wave state and magnetic transition Mi contribution only. The 
M\ amplitude is calculated up to N 2 LO. Three-nucleon forces are needed up to N 2 LO order for cut- 
off independent results. The triton binding energy and nd scattering length in the triton channel 
have been used to fix them. Hence the The photon polarization parameter in total is determined as 
R c = -[0.387(LO) + 0.016(iVLO) + 0.009(N 2 LO)] = [-0.412 ± 0.003]. This converges order by order 
in low energy expansion and also is cut-off independent at this order. We notice that our calculation 
has a systematic error which is now smaller than the experimental error bar. 

VI. ACKNOWLEDGMENTS 

I would like to thank Harald W. Griefihammer for useful and valuable comments and L. Marcucci for 
providing the recent theoretical results. This work was performed under the auspices of the university 



8 



of Arak under contract No. 11-18680. 



[1] M.W. Konijnenberg et al., Phys. Lett. B 205, 215 (1988). 

[2] J. Golak et al, Phys. Rev. C 62, 054005 (2000). 

[3] R. Skibinski et al, Phys. Rev. C 67, 054001 (2003). 

[4] W. Schadow, O. Nohadani, and W. Sandhas, Phys. Rev. C 63, 044006 (2001). 

[5] A. Deltuva, L. P. Yuan, J. Adam, Jr., A. C. Fonseca, and P. U. Sauer, Phys. Rev. C 69, 034004 (2004). 

[6] J. Carlson and R. Schiavilla, Rev. Mod. Phys. 70, 743 (1998). 

[7] M. Viviani, R. Schiavilla, and A. Kievsky, Phys. Rev. C 54, 534 (1996). 

[8] L.E. Marcucci, M. Viviani, R. Schiavilla, A. Kievsky and S. Rosati, Phys. Rev. C 72, 014001 (2005). 

[9] H. Sadeghi and S. Bayegan, Nucl. Phys. A 753, 291 (2005). 

[10] H. Sadeghi and S. Bayegan, H. W. Griefihammer, Phys. Lett. B 643, 263 (2006). 

[11] E.T. Jurney and H.T. Motz, Argonne National Laboratory Report ANL-6797, 1963, P. 236 (unpublished). 

[12] J.S. Merritt, J.G.V. Taylor and A.W. Boyd, Nucl. Sci. Eng 34, 195 (1968). 

[13] E.T. Jurney, P.J. Bendt and J.C. Browne, Phys. Rev. C 25, 2810 (1982). 

[14] P. F. Bedaque, G. Rupak, H. W. Griefihammer and H.-W. Hammer, Nucl. Phys. A 714, 589 (2003). 

[15] H. W. Griefihammer, Nucl. Phys. A 744, 192 (2004). 

[16] J.-W. Chen, G. Rupak, M.J. Savage, Phys. Lett. B 464, 1 (1999). 

[17] V. Bernard, T. R. Hemmert and U. G. Meifiner, Nucl. Phys. A 732,149 (2004). 

[18] E. Epelbaum, W. Gloeckle and U. G. Meifiner, Eur. Phys. J. A 19, 125 (2004). 



9 



